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1  .   Introduction 
It  has  been  common  practice  in  the  recent  Real  Business  Cycle  (RBC) 
literature  to  approximate,  quadratical ly ,  the  return  function  about  the 
steady  state  and  then  to  use  this  approximate  return  function  as  the 
basis  for  generating  the  economy's  equilibrium  time  series.   This  is 
done  for  well  known  reasons  of  analytic  and  computational  simplicity: 
with  a  quadratic  return  function  the  decision  rules  are  linear  and  may 
be  easily  determined.   To  compute  the  optimal  decision  rules  numerically 
via  standard  value  iteration  procedures  is  simply  too  intensive  when  the 
number  of  decision  and  state  variables  is  large.   Nevertheless,  it  is 
legitimate  to  question  the  extent  to  which  accuracy  is  compromised  using 
such  approximate  procedures. 

The  answer  to  this  question  may  very  well  be  entirely  model  spe- 
cific.  Indeed,  if  it  were  feasible  to  undertake  such  a  comparative 
evaluation  for  all  model  contexts,  there  would  be  no  need  to  resort  to 
such  approximations  in  the  first  place.   In  this  paper  we  undertake  such 
a  comparison  for  the  basic  one  good  stochastic  growth  model  as  analyzed 
in  Hansen  (1985)  (his  'divisible  labor  economy').   That  is,  we  solve  for 
the  optimal  decision  rules  and  equilibrium  time  series  using  both  stan- 
dard value  iterative  and  quadratic  approximation  techniques  and  compare 
results.   Some  familiarity  with  the  basic  methodology  is  assumed. 

2.   The  Basic  Model 
The  model  we  investigate  is  a  familiar  one 


oo 

,t 
ma 


x  E(  I  pu  u(c   1-N  )) 
t=0 


s.t.  Ct  +  Zt  -  '^t^t^t'  k<)  Siven 

k    =  zt  +  (1  -  fi)k   ,  0  <  N   <  1  ,  where 

c  ,  k  ,  N   and  z   arc,  respectively,  per  capita  consumption,  capital 
stock,  labor  service  supplied,  and  investment  in  period  t,  E  is  the 
expectations  operator,  p  the  period  discount  factor,  Q  the  depreciation 
factor,  U(  )  the  representative  agent's  period  utility  function,  f(  ) 
the  production  function,  and  A   the  random  period  shock  to  technology. 
The  precise  specification  of  the  functional  forms  employed  here  is  as 

follows:   U(ct,l-Nt)  =  glc^l-N^1"*]6,  and  f(kt,Nt)  =  \    ["n*~V\. 
For  v  =  L  =  1  and  a   =  .36,  the  above  technology  specification  is  the 
same  as  Hansen  (1985);  the  preference  specification  is  taken  from  Kydland 
and  Prescott  (1982)  . 

The  choice  of  shock  process  requires  a  somewhat  more  extensive 
discussion.   In  Hansen's  (1985)  work,  the  shock  to  technology  is 
described  by  first  order  Markov  process  of  the  form 

(2)  .  At+l  =  PAt  +  et 

with  p  =  .95  and  e   lognormally  distributed  with  Ee   =  1  -  p  =  . 05  and 

a  =  .00712.   While  such  a  shock  process  is  perfectly  natural  for  the 

quadratic  return,  linear  decision  rule  setting,  it  does  not  immediately 

translate  into  the  standard  value  iterative  context.   Under  these 

latter  procedures,  it  is  first  necessary  to  define,  precisely,  the 

region  in  which  the  state  variables  --  in  this  case,  capital  stock  and 

the  shock  to  technology  --  lie  over  the  problem's  infinite  horizon.   For 

Hansen's  specification  of  the  production  technology,  the  range  of 

possible  capital  stock  values  will  be  bounded  above  by  some  k    provided 
r  J  max 

the  shock  to  technology  is  itself  bounded. 


Given  upper  and  lowci  bounds  on  the  shock,  it  is  then  necessary  (sec, 
e.g.,  DanLhine  ami  Donaldson  (1981),  Mchra  and  Prescott  (1985),  and 
Greenwood  and  Huffman  (1986))  to  select  a  discrete  set  of  N  possible 
values  of  A,  typically  with  mean  one,  the  probability  ol  relative  occur- 
rence of  which  is  governed  by  a  prespecified  N  dimensional  probability 
transition  matrix. 

Hansen's  (1985)  shock  process,  however,  is  not  uniformly  bounded 
above  over  the  infinite  horizon  and  thus  has  no  direct  representation  in 
the  transition  matrix  setting.   In  order  to  undertake  our  comparisons, 
we  therefore  specify  a  two  state  Markov  process  on  the  shock  to  tech- 
nology which  has  the  same  mean  and  variance  as  the  process  chosen  by 
Hansen  (1985).    Thus,  for  both  the  value  iterative  and  quadratic 
approximate  procedures  we  assume  that  the  Markov  process  on  the  tech- 
nology shock  assumes  values  A  =  .98  and  A  =  1.02  with  transition 
density 


K         =  .98      1.02 


A  =   .98 
1.02 


.97 
.03 


.03 
.97 


We  denote  the  conditional  distribution  of  next  period's  shock  given 
today's  shock  of  A  by  dF(A  ;A) . 


3.   The  Standard  Solution 
The  sequence  of  approximating  value  functions  is  described  by  the 
recursive  equation 


(3) 

2 

V  (k,A)  =  max{U(f(k,N  )A  -  k  +  (l-tt)k,l-N  )  +  p  1      V   ,(k  ,A  )d¥(\    ;\)} , 
n  n      n  n  n-lntt 

s.t .  0  <  N    <  ] 

n 

(l-fi)k  <  k   <.  f(k,N  )A  +  (l-fi)k 
n  -      n 

where  the  subscript  n  denotes  the  nth  iteration.   Experience  has  shown  us 

that  optimizing  over  next  period's  capital,  rather  than  investment, 

2 
substantially  reduces  round  off  errors. 

The  state  and  decision  variables  are  constrained  to  assume  values 
in  a  discrete  set  or  "grid"  which  reasonably  approximates  their  actual 
domain  of  definition.   A  simple  serial  search  procedure  for  determining 
the  optimal  (k  ,N  )  by  evaluating  the  above  expression  for  all  possible 
(k  ,N  )  is,  however,  computationally  very  slow.   Assuming  as  many  as  250 
distinct  possible  values  for  the  capital  stock,  two  values  of  the  tech- 
nology shock  and  100  possible  levels  of  labor  supply,  each  stage  of  the 
iterature  process  would  require  (250  x  2)  x  (250  x  100)  (number  of 
possible  state  variable  combinations  x  number  of  possible  pairs  of 
decision  levels)  =  12,500,000  independent  evaluations  of  the  above 
expression  to  determine  the  maximum.   Although  this  number  can  be  sub- 
stantially reduced  by  applying  known  properties  of  the  decision  variables 
(such  as  the  fact  that  k   is  increasing  in  k) ,  in  light  of  the  hundred 
plus  iterations  that  may  be  necessary  for  policy  convergence,  the  total 
computation  time  is  f rustratingly  large. 

A  natural  device  for  reducing  these  computations  by  a  factor  of  100 

is  to  devise  a  simple  and  rapid  procedure  for  determining  directly  the 

optimal  N   associated  with  each  possible  choice  of  k  .   This  is  what  we 
n  n 

do.   Noting  that,  for  a  given  (k,X): 


U)  searching  over  the  region  of  ordered  pairs  (k  ,N  )  for  which 

n   n 

0  ^  N   ^  1  and  (1  -  Q)k  v  k   v  Hk.N  )  +  ( 1  -  fi)k  is  equivalent  to 
n  n  -      n  M 

searching  ovrr  the  region  (1  -  Q)k  <  k   <  f(k,l  )  +  (1  -  Q)k  and 

n  -      n 

] 


N0  <  N   <  1  where  N0  =  N0(k  ,k,A)  s 


[v|kn  -  (]  -  Q)k]"|(l-a)v 
1       ,  ant 


(ii)  maximizing  U(  )  +  V    (  )  with  respect  to  N   is  equivalent 

to  maximizing  U(  )  alone  with  respect  to  N  ,  since  V    (  )  is 

independent  of  N  , 
n 

we  can  solve  for  the  optimal  N  associated  with  each  (k  ,k,A)  triple  by 
maximizing 


(4)  (aN(1-a)v  -  b)*6(l  -  N  )(1"*)6c 

n  n 


LkaVA 

with  respect  to  N  ,  where  a  =  a(k,\)  =  >  0,  b  =  b(k  ,k)  = 

r        n  '       v   —  '        n 

k  -  (1  -  Q)k  >  0,  and  c  =  t.   Since  equation  (A)  is  a  concave  function 

of  N  we  may  differentiate  and  set  the  resultant  expression  equal  to  zero 

to  determine  the  optimal  N  .   This  yields  the  expression 


aNv    '      -  b  1  -  N 

n  n 


Rearranging  terms  this  may  be  rewritten  as 


(6)  AN(1"a)V  -  BN*1-0^-1  -  C  =  0 

n         n 


where  A  =  (l-v)a  +  va(l-a)u  >  0,  B  =  ya(l-a)v  >  0  and  C  =  (l-y)b.   Rear- 
ranging equation  (6)  again  gives 


l-(l-a)\> 


(7) 


AN   -  B  -  CN 
n  n 


0,  or 


(8) 


n    A    An 


Equation  (8)  is  amenable  to  standard  fixed  point  procedures 
Graphing  each  side  of  (8),  as  per  Figure  (1)  below,  it  is 


B 
A 


B  Jk  CMl-(l-a)v 
A   An 


Figure  (1) 


app 


arent  that  the  recursive  iteration  procedure  defined  by 


N(0)  E  1 
n 


nO)  =  B  +  C  (0)  l-(l-a)v 
n     A   A  n 


(9) 


,(2)  =  B  +  C   (1)  l-(l-a)v 
n     A   A  n 


N(k)  =  B  +  C  (k-D.l-Cl-a^ 

n     A   A   n 


will  approximate  the  solution.   The  remaining  issue  is  to  establish  a 
convergence  criterion  for  the  procedure.   Suppose  the  solution  to 

equation  (8)  is  an  N   =  jj  +  [   N1-(l-a)v   thcn  f   gn  N0) 

n    A    A   n  n 


(10) 


N(i+1)  =  B  +  C   (i)  l-(l-a)v  =  B  +  C    l-(l-a)v 
n       A   A   n  A   A   n 


+  ^(l-(l-a)v)N-(1-a)V(N(i)  -  N  ) 
A  n        n      n 


=  N  +  £(l-(l-a)V)N_(1"o)v(N(i)  -  N  ) 
n   A  n        n      n 


If  we  desire  a  solution  for  which  N     -  N   <  .5||K||,  where  IlKll  denotes 

n      n 

the  width  of  the  partition  of  possible  capital  stock  levels,  then  the 
convergence  requirement  becomes: 


H(i+D  .  N(i) 
n        n 


N  +  ^(l-(l-a)v)N-(1-a)V(N(i)  -  N  )  -  N(i) 
n   A  n        n      n     n 


[^l-(l-a)v)N-(1"a)V  -  l](N(i)  -  N  ) 
A  n  n      n 


<   T(l-(l-a)v)N       -  1 
-   A  n 


(.5  ||K||),  where 


-7,  and  N   can  be  estimated  from  the  steady  state  certainty  levels  of 

An  J  J 

■      ,    >  ,  x.     v    .   C  ~  (1  -  v)[kss  -  (1  -  fi)kss]v  „  ~  .   „ss 
capital  and  labor;  that  is  -  = u-i — ,  N   =  .3  =  N   , 

[(1  -  Y)  +  Yd  "  a)v]Lk  A. 
and  (1  -  (1  -  y)u)  =  -3°  (in  the  case  of  Hansen's  (1985)  choice  of  param- 


eters).  Given  the  determination  of  the  optimal  N   =  N  (k  ,k,X),  the 

related  optimal  k  was  determined  via  a  modified  exhaustive  search 
n 

procedure.   The  time  series  output  of  the  model  was  constructed  by  first 
generating  a  sequence  of  three  thousand  shocks  to  the  technology  with 


respect  to  the  chosen  transition  matrix  and  then  by  allowing  the  economy 
to  evolve  from  an  arbitrarily  determined  initial  state  via  repeated 
application  of  the  optimal  policy  rules.   We  have  not  subjected  the  data 
to  any  detrending  procedure  as  we  thought  it  worthwhile  first  to  assess 

the  model's  direct  output  without  the  modifications  a  smoothing  pro- 

3 
cedure  would  introduce.    Limiting  ourselves  for  the  moment  to  the  sta- 
tistics reported  in  Hansen  (1985),  we  find  the  time  series  properties  of 
this  optimal  economy  to  be  summarized  below: 


Table  (1) 

B  =  .96,  a   =  .36,  y  =  -33,  ft  =  .1,  \j  =  1.02,  A2  =  .98, 
Persistence  =  .97,  L  =  1 
Capital  Stock  Partition  =  .005 


(a)  standard  deviation     (b)  correlation  with  output 
(percent) 

output:  2.43  1.00 

consumption:  2.09                         .898 

investment:  4.49                         .881 

capital  stock:  2.22                        .868 

hours:  .69                        .99 

productivity:  2.23                        .276 


We  note  that  for  a  state  space  of  250  possible  capital  stock  levels, 
and  100  possible  labor  supply  levels,  the  entire  computational  procedure 
required  31.33  minutes  of  VAX-11-780  CPU  time. 

3.   An  Approximate  Solution 
Our  procedure  is  similar  to  that  of  Section  (2)  except  that  we 
undertake  a  quadratic  approximation  of  the  return  function  about  the 
economy's  certainty  steady  state.   As  before,  we  first  outline  our 
iterative  procedure  and  then  present  our  comparative  time  series  results. 


Let  the  superscript  ss  denote  tlir  steady  state  certainty  (A    1) 
value  of  the  relevant  variable.    Following  Kydland  and  Prescott  (  1982) 
and  Hansen  (198ri),  we  approximate  the  agent's  utility  function  about  the 
steady  state  by  the  following  expression: 


(11) 


u"(k,N,A,z)  =  U  +  BX  +  X'QX,  where 


(12) 


X  = 


*2 

*3 

x4 


-  N 


k 
N 
X 
z  -  z 


ss 


ss 


ss 


ss 


(13)  B  = 


I*  /.  SS      .    -_SS  •  SS    SS  ■.      a  ,  SS     ,    ..SS  -.SS    SS'.  i   1 
u(k   +  Ak,N   ,A   ,z   )  -  u(k   -  AR,N   ,A   ,z   )]^- 

k 


■  a  .  SS  ..SS    .   .  SS   SSN    a  ,  ss  ..ss    .   .ss   ss .  .  1 
u(k   ,N   +  An,A   ,z   )  -  u(k   ,N   -  A^A   ,z   )]^- 


a,,  SS  .-SS  .SS     A     S  S       a  /,  SS  «TSS  \  SS      A     Sb\  i   1 

u(k   ,N   ,A   +  AA,z   )  -  u(k   ,N   ,A   -  A^,z   )]^~ 


a,.SS  ..ss  ,SS   SS     .  N     a  .  ss  ..ss  .ss   SS     .  n.  i  1 

u(k   ,N   ,A   ,z   +  A  )  -  u(k   ,N   ,A   ,z   -  A,)]^" 


where  G(k,N,A,z)  =  u(f(k,N)A  -  z,  1  -  N) ; 


(14) 


U=u(k   ,N   ,A   ,z   );  and 
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y   is    a    4    y    A   matrix   with    diagonal    entries    oi    the    form: 


qu    =    |u(k        +Ak'N      <A      »z      )    -    u(k      ,N      ,A      ,z      ) 


'  /I   SS  A  »TSS      iSS  ss,       1 

+   u(k        -    AN      ,A      ,z      ) — r 

k  2A: 

k 


(15) 


I^/.SS    .TSS     *  SS       SS  .     >.  a  .    ss    ..ss     ,  s  s        ss  >. 

u(k      ,N      ,A      ,z        +  A   )    -   u(k      ,N      ,K      ,z      ) 


+   -   G(kSS,NSS,ASS,zSS    -   A   )]~ 

Z      2A^ 
k 


and   off   diagonal   elements   of   the   form,    e.g. 


q13  =    [G(kSS   +  Ak,NSS,ASS   +  AA,zss)    -   G(kSS    -   Ak,NSS,ASS   +  AA,zSS) 


(16) 


A/i  SS  A         .TSS     .SS  .  SSv  m  #,  SS  A         ,,SS      .SS  .  S  S  ■   -,  1 

u(k        +  A    ,N      ,\       -   A    ,z      )    +  u(k        -  A,  ,N      ,A        -   A.,z      )] 


k,.,      ,»  ^A,^      ,_      uv.  „k,„      ,„  „A,«.      „  •   A.)8 

k  A 


In  the  above  expressions,  the  symbol  A  indicates  the  deviation  from  the 
steady  state  while  U  represents  the  certainty  steady  state  level  of 
period  utility.   The  deviations  themselves  were  chosen  proportional  to 
the  steady  state  levels  with  the  same  proportionately  constant  across 
all  state  variables. 

Turning  to  the  iterative  procedure  itself,  it  is  evident  that  for 
fixed  (k,A),  V^(k,A)  can  be  expressed  as 

(17)     vj(k,\)  =  CONSTANT  +  b^  +   I  (q   +  q.2)x.x2  +  q^x* 

where  the  constant  has  the  value 


n 


(18)      CONSTANT  -  u  +  1   h   x  .    ■*      lq..x.x..  and  X.  =  (1  -  U)V,    -   zBS 

i*2  '  '    i*2  1-1  3  J       U 
j*2 


The  superscript  A  is  intended  to  remind  us  that  these  are  approximate 

value  functions.   We  note  also  that  the  precise  represent  at  ions  of  the 

x.'s,  b.'s  and  q.'s  have  been  suppressed  for  ease  of  presentation. 
i      l        ij 

Noting  that  equation  (17)  is  of  the  form  V  =  C  +  Bx2  +  Ax£,  it  is 
clear  that  the  maximum  will  occur  at 


~     s  s 
from  which  the  optimal  N0  =  N   +  x   is  easily  computed.   It  follows  that 

(20)  v{(k,A)  =  u(f(k,N0)A,l  -  N0)  for  all  (k,A). 

More   generally, 


(21)      Ak,A)   =.         max  {V*(k,N   ,A,k   )   +   P  I  vj      (k     A   )dF(A     A)} 

0<N<1  t=l 


(1    -  fi)k  <   k     <   f(k,N   )A  -    (1  -  fi)k 
-  .  n  -  n 

We  point  out  that  in  this  expression  we  optimize  over  next  period's  capi- 
tal stock,  rather  than  this  period's  investment.   This  requires  some  mod- 
ification of  U  (  ) ,  but  again  seems  to  substantially  reduce  errors 
resulting  from  rounding  off. 

An  examination  of  the  region  defining  the  set  of  possible  choices  of 
the  state  variables  allows  us  to  express  (21)  equivalently  as 
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A  *  2   A 

(22)  V  (k,A)  =    max     u  (k,N  ,A,k  )  *  0  I   V   n(k  ,A  )dF(A  ,A)} 

n  n     n  n-l   n   t      t 

s.t  .     (1  -  fi)k  <  k   <  f Ck.l)  +'  (1  -  fi)k 

n  - 

N„(k  ,k,A)  <  N   <  1 
n      -  n  - 

Substituting  and  rearranging  terms  gives: 

(23)  v\k,A)  =   max  {     max  {C(k,A,k  ) 

(l-fi)k  <  k  <  f(k,N  )  +  (l-O)k   N0(k  ,k,A)  <  N   <  1        n 
—  n  -      n  n      —  n  — 

+  V2  +  .J2(qj2  +  q2j)XjX2  +  q22X2}}  Where 

(24)  C(k,A,k  )  =  u(kSS,NSS,ASS,zSS)  +  I   b.x.  +  I   q..x.x. 

it2   ii    lt2    U  i  J 

j*2 


+  M  v*  (k  ,\  )dF(\  ,\) 

t  =  l 


As  before  the  optimum  N  can  be  computed  directly  by  differentiating 
the  expression  of  line  (24)  to  give 


-0>2  +   1   2q  x.) 

*  -  -  £ Jt^ '  R°" NSS 


We  note  also  that  if  N   >  1,  then  N  =  1;  if  N   <  N0 ,  N   =  N0  and  if 

n  n  n        n 

No  <  N   <  1,  then  N   =  N  .   Given  this  information,  the  computation  of 
n  n    n 

the  optimal  k   reduces  to  solving 


(26)   V^k.A)  =    max    {u"(k,N  ,X,k  )  +  p  I  V^   (k  ,A  )dF(A  ,A)  ] 
(1  -  D)k  <  k  <  f(k,l)  +  (1  -  fi)k 
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This  latter  maximization  was  solved  using  a  modi  lied  standard  grid  search 
procedures.   The  resulting  savings  policy  function  is  linear  and  is  the 
same  as  li  this  expression  had  also  been  solved  directly.   Although  we 
could  as  well  directly  solve  for  the  capital  next  period,  the  applica- 
tion of  the  linear  decision  rule  so  obtained  would  not  necessarily  con- 
fine the  economy  to  the  prechosen  state  space.   In  effect,  the  discretiz- 
ing  of  the  state  space  precludes  us  from  using  the  Kydland  and  Prescott 
(1982)  procedure.   The  results  of  the  simulation  for  the  same  parameters 
as  in  Table  (1)  are  reported  below: 


Table  (2) 

p  =  .96,  a  -    .36,  v  =  -33,  fi  =  .1,  \t  =   1.02,  A2  =  -98, 

Persistence  =  .97,  A  =  .001 

Capital  Stock  Partition  =  .005 


(a)  standard  deviation     (b)  correlation  with  output 
(percent) 

output:                    2.18  1.00 

consumption:                2.27  .992 

investment:                2.6A  .678 

capital  stock:               .59  .899 

hours:                      .31  .99 

productivity:               2.02  .062 


Comparing  Tables  (1)  and  (2)  we  note  that  the  standard  deviations  of 
the  output,  consumption,  hours  and  productivity  series  are  more  or  less 
the  same.   The  exact  solution  has  substantially  greater  standard  devia- 
tion of  investment  and  capital  stock,  however.   It  should  be  noted  that 
each  of  the  time  series  from  which  these  statistics  were  calculated  was 
(1)  subjected  to  identical  shocks  (randomly  generated)  and  (2)  had  its 
origin  in  the  same  initial  capital  stock  level.   It  is  also  worth  noting 
that  the  levels  of  all  series  (expected)  were  nearly  identical. 
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4 .      Sensitivity  Analysis 

There  are  a  number  of  related  issues  which  ought  to  be  addressed  if 
we  are  to  assess  fully  the  appropri ateness  of  this  approximation.   In 
particular,  we  examine  (1)  the  accuracy  of  the  approximation  with  regard 
to  the  magnitude  of  the  deviation  from  steady  state  levels,  (2)  the 
comparative  time  paths  to  the  steady  state  of  the  approximate  and  exact 
economies  and  (3)  the  sensitivity  of  the  approximation  to  changes  in 
parameter  values.   Each  is  considered  in  turn: 

4.1   Magnitude  of  Deviations:   As  noted  in  Section  4,  for  the 
approximate  results  reported  thus  far,  all  deviation  from  the  steady 
state  (the  A's)  were  chosen  proportional  to  the  corresponding  steady 
state  level  itself  with  the  proportionality  factor  set  at  .001.   In  this 
subsection  we  recompute  the  statistics  of  Table  (2)  for  the  case  in  which 
this  proportionality  factor  is  .01  and  .02.   These  results  are  reported 
below: 


Table  (3) 

A  =  .01 
a  =    .36,  p  =  .96,  v  =  .33,  fi  =  .1,  \a  =  1.02,  X2   =    -9i 
Persistence  =  .97,  L  =  1 
Capital  Stock  Partition  =  .005 


(a)  standard  deviations  (b)  correlation  with  output 
(percent) 

output:                    2.87  1.00 

consumption:                2.51  .968 

investment:                 7.43  .879 

capital  stock:              3.36  .895 

hours:                    1.12  .986 

productivity:               2.43  .350 

A  =  .02 
(other  parameters  unchanged) 
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(a)  standard  deviations    (b)  correlations  with  output 

output :  3.40  ] .00 

ronsumpt 1  on :  3.05  .956 

investment :  9.16  .898 

capatal  stork:  5. A]  .910 

hours  :  ]  .  33  .  9(>n 

productivity:  2.87  .523 


These  tables  reveal  the  tendency  of  the  standard  deviation  in  all 
series  to  increase  as  the  proportional  deviation  A  increases;  capital 
stock  variation  appears  to  increase  at  the  greatest  rate.   This  is  not 
surprising  as  the  quadratic  approximation  is  effectively  reducing  the 
level  of  concavity  in  the  utility  function  which  has  the  consequence  of 
making  increased  variation  in  all  the  series  optimal.   With  regard  to  the 
investment  function  in  particular,  its  linearity  under  the  quadratic 
approximate  regime  contrasts  with  its  concavity  under  the  'exact'  formu- 
lation.  For  any  given  variation  in  output,  greater  variation  in  invest- 
ment will  be  observed  under  the  linear  decision  rule. 

We  also  note  that  the  correlation  of  both  consumption  and  hours  with 
output  declines  slightly  with  the  increased  A,  while  wage  (productivity) 
correlation  with  output  increases  very  dramatically.   The  correlation  of 
output  with  the  remaining  other  series  is  substantially  unaffected.   We 
are  unable  to  provide  intuition  as  to  why  this  should  be  so. 

A. 2   Convergence  to  Steady  State:   As  noted  earlier,  our  program  is 
designed  so  that  all  time  series  are  generated  by  subjecting  the  tech- 
nology to  the  exact  same  sequence  of  random  shocks.   The  initial  capital 
shock  level  is  identical  across  all  cases  as  well.   This  set-up  has  the 
advantage  of  allowing  us  to  study  relative  convergence  rates  under  both 
procedures . 
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In  general,  the  exact  economy  converges  to  the  steady  state  capital 
stock  level  much  less  rapidly  than  does  the  approximate  economy.   The 
exact  behavior  is  as  follows:   both  economies  linger  at  the  boundary  of 
the  steady  state  without  actually  entering  it  for  many  periods.   This 
may  simply  reflect  a  long  sequence  of  the  poor  random  shocks.   However, 
the  approximate  economy  reaches  the  boundary  in  approximately  one-fourth 
the  time  required  for  the  exact  economy.   This  is  all  the  more  striking 
when  we  realize  that  the  range  of  the  stationary  distribution  on  capital 
stock  for  the  exact  economy  fully  contains  the  range  for  the  approximate 
economy.   These  figures  are,  respectively,  [1.2,  1.26]  and  [1.22,  1.235]. 
Thus,  not  only  does  the  exact  economy  require  more  time  to  converge,  but 
also  the  amount  by  which  capital  must  increase  to  enter  the  steady  state 
(ko  <  1.2)  is  less.   These  reported  range  values  are  also  consistent  with 
the  greater  relative  variability  of  capital  stock  for  the  exact  economy. 

4.3   Sensitivity  to  Parameters:   Of  all  the  parameters  necessary  to 
define  this  class  of  problems,  the  choice  of  the  preference  parameter  y 
is  perhaps  most  open  to  question  and  we  examine  the  sensitivity  of  our 
approximation  to  changes  in  that  value.   This  has  the  added  advantage  of 
leaving  the  steady  state  basis  of  the  approximation  unchanged.   We  find 
that  with  the  A  =  .01  deviation,  the  discrepancies  between  the  exact  and 
approximate  cases  can  vary  substantially  with  the  choice  of  preference 
parameters.   This  is  illustrated  below  where  we  first  allow  for  changes 
in  v  (Table  4) . 
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A,  - 


Tabic  4 

B  =  .96,  6=  -.5,  fi  =  .1,  A  =  .01 

1.02,  AL.  =  .98,  Persistence  =  .97,  L 
Capita]  Stock  Partition  =  .005 
Y  =  -25 


output : 
consumption: 
investment : 
capital  stock: 
hours : 
productivity: 


Exact  Economy 
(a)  standard  deviation 

(a)  (b) 


2.56 
2.12 
6.32 
2.49 
1.13 
2.33 


1.00 
.971 
.875 
.857 
.991 
.286 


Approximate  Economy 
(b)  correlation  with  output 


(a) 

3.63 
4.16 
14.52 
6.23 
1.64 
3.16 


(b) 


.908 
.823 
.931 
.971 

.597 


output: 
consumption: 
investment : 
capital  stock: 
hours : 
productivity: 


Y  =  .40 

Exact  Economy 
(a)  (b) 


2.48 
1.99 
5.80 
2.47 
.82 
2.28 


1.00 
.976 
.888 
.870 
.994 
.307 


Approximate  Economy 
(a)  (b) 


2.89 
2.36 
5.54 
2.89 
.96 
2.32 


1.00 
.981 
.909 
.879 
.988 
.297 


5 .   Concluding  Comments 

The  figures  reported  here  suggest  that  the  quadratic  approximation 
will  provide  a  very  reasonable  approximation  to  its  exact  counterpart 
when  the  shocks  are  low  and  the  deviations  which  form  the  basis  of  the 
approximation  are  small.   It  is  in  the  capital  stock  and  investment 
series  that  the  greatest  discrepancy  arises. 

All  of  this  is  a  bit  unsatisfying,  however,  as  there  are  many 
sources  of  approximation  error  in  a  numerical  exercise  as  complex  as 
this  and  it  seems  impossible  to  assign  any  discrepancies  in  the  'exact' 
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and  'approximate'  solutions  to  a  particular  source.   Thus  the  use  of 
such  methods  in  increasingly  complex  numerical  exercises  will  remain, 
invariahly,  partially  an  act  of  faith. 


19 


FOOTNOTES 

1.  This  was  accomplished  using  the  method  of  moments. 

2.  If  we  optimize  directly  over  next  period's  capita]  stock  then 
investment  z  ami  consumption  are  computed  as  residuals.   If  we 
optimize  over  z,  the  z  is  discretized,  (,l-fl)k  is  discretized  and 
the  sum  z  +  (l-ft)k  discretized  --  which  gives  three  sources  of 
round  off  errors  in  the  measurement  of  the  critical  next  period's 
capital  stock  level. 

3.  Since  the  output  data  of  this  model  does  not,  by  construction, 
exhibit  any  trend,  it  is  not  entirely  clear  that  such  a  procedure 
would  be  appropriate  anyway. 

4.  These  steady  state  values  were  computed  as  in  Kydland  and  Prescott 
(1982). 
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